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ABSTRACT 

We investigate the effect of primordial non-Gaussianity of the local /nl type on 
the auto- and cross-power spectrum of dark matter haloes using simulations of the 
ACDM cosmology. We perform a series of large N-body simulations of both positive 
and negative /nl, spanning the range between 10 and 100. Theoretical models predict 
a scale-dependent bias correction A6(fc,/NL) that depends on the linear halo bias 
b{M). We measure the power spectra for a range of halo mass and redshifts covering 
the relevant range of existing galaxy and quasar populations. We show that auto and 
cross-correlation analyses of bias are consistent with each other. We find that for low 
wavenumbers with k < 0.03 hMpc~^ the theory and the simulations agree well with 
each other for biased haloes with b{M) > 1.5. We show that a scale-independent bias 
correction improves the comparison between theory and simulations on smaller scales, 
where the scale-dependent effect rapidly becomes negligible. The current limits on /nl 
from Slosar et al. (2008) come mostly from very large scales k < 0.01 hMpc~^ and, 
therefore, remain valid. For the halo samples with b{M) < 1.5 — 2 we find that the scale- 
dependent bias from non-Gaussianity actually exceeds the theoretical predictions. Our 
results are consistent with the bias correction scaling linearly with /nl- 
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1 INTRODUCTION 



, Generic inflationary models based on the slow roll of a scalar 
5_j ' field predict a nearly scale-invariant and Gaussian spectrum 

I of primordial curvature fluctuations (see Bartolo et al. 2004 
for a review). While the latest measurements of the cos- 
mic microwave background (CMB) anisotropies favour a 
slightly red power spectrum (Komatsu et al. 2008), no sig- 
nificant detection of primordial non-Gaussianity has been 
reported as yet from CMB and large-scale structures mea- 
surements. Nevertheless, improving the current limits would 
still strongly constrain mechanisms for the generation of cos- 
mological perturbations. 

Non-Gaussianity can be generated by nonlinearities in 
the relation between the primordial curvature perturbation 
and the inflaton field (e.g. Salopek & Bond 1990; Gangui et 
al. 1994), interaction of scalar fields (e.g. Falk et al. 1993) or 
deviation from the (Bunch-Davies) ground state (e.g. Les- 
gourgues et al. 1997). A wide class of inflationary scenarios 
lead to non-Gaussianity of the local type, which depends on 
the local value of the potential only. In these models, de- 
viation from Gaussianity can be conveniently parametrised 



by a noiflinear coupling parameter /nl through the relation 
(e.g. Komatsu & Spergel 2001) 



- /nl {H^f - {4>(xf)) 



(1) 



where 0(x) is the Gaussian part of the curvature perturba- 
tion in the matter area. While single inflaton scenarios pre- 
dict /nl much less than unity, multi-fleld inflation models 
can generate /nl 2> 1 (Linde & Mukhanov 1997; Lyth et al. 
2003; CreminelU 2003; Dvah et al. 2004; Zaldarriaga 2004; 
Arkani-Hamed et al. 2004; Alishahiha et al. 2004). Alterna- 
tives to inflation, such as cyclic /ekpy rot ic model also predict 
large non-Gaussianity of local type (Creminelli & Senatore 
2007, Buchbinder et al. 2008, Lehners & Steinhardt 2008). 

Higher order statistics of the curvature perturbation 
such as the bispectrum can be computed straightfor- 
wardly from a perturbative expansion of the homogeneous 
Robertson- Walker background (e.g. Acquaviva et al. 2003; 
Maldacena 2003). These statistics are related to those of the 
CMB temperature anisotropy through the radiation trans- 
fer function, which can be computed accurately using, e.g., 
CMBFAST (Seljak & Zaldarriaga 1996). Thus far, analysis of 
the CMB bispectrum indicates that the data are fully con- 
sistent with Gaussianity, with |/nl| <,W0 (Komatsu et al. 
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2003; Creminelli et al. 2007; Komatsu et al. 2008; Smith et 
al. 2009; see, however, Yadav & Wandelt 2008 who report 
a detection at the 2.5(t level), providing strong evidence for 
the quantum origin of the primordial fluctuations. 

Large-scale structures offer another route to test for 
the presence of primordial non-Gaussianity. It has long been 
recognised that departure from Gaussiamity can significantly 
affect the high mass tail of the dark matter halo distribution 
(Lucchin & Matarrese 1988; Colafrancesco et al. 1989; Chiu 
et al. 1998; Robinson & Baker 2000; Matarrese et al. 2000; 
Mathis et al. 2004, Kang et al. 2007; Grossi et al. 2007). Fol- 
lowing this approach. X-ray cluster counts have been used 
to constrain the amount of non-Gaussianity (e.g. Koyarna 
et al. 1999, Robinson et al. 2000; Willick 2000; Amara & 
Refregier 2004). Galaxy clustering is also sensitive to the 
statistical properties of the primeval fluctuations. Indeed, 
Grinstein & Wise (1986) pointed out early that primordial 
non-Gaussianity could significantly increase the amplitude 
of the two-point correlation of galaxies and clusters on large 
scales. However, recent work has mostly focused on higher 
order statistics such as the bispectrum (Scoccimarro et al. 
2004; Sefusatti & Komatsu 2007). 

Dalai et al. (2008) have recently sparked renewed in- 
terest in the clustering of rare objects by demonstrating 
the strong scale-dependent bias arising from primordial non- 
Gaussianity of the local type. It can be shown that the latter 
contributes a scale-dependent bias of the form (Dalai et al. 
2008; Matarrese & Verde 2008; Slosar et al. 2008) 

A64fc, /nl) = [6(M) - 1] 5c^5^^|^ , (2) 

where 6(M) is the linear bias parameter, Hq is the Hubble 
parameter, T{k) is the matter transfer function, D{z) is the 
growth factor normalised to (1 -I- z)~^ in the matter era and 
5c ^ 1.68 is the present-day (linear) critical density thresh- 
old. While the derivation of this non-Gaussian bias correc- 
tion presented in Dalai et al. (2008) and Matarrese & Verde 
(2008) is strictly valid only for the highest peaks of the den- 
sity field, the peak-background split argument invoked by 
Slosar et al. (2008) suggests that cq. (2) should apply to all 
peaks unrestrictedly, but is only valid in the limit of long 
wavelength modes so that the background can be approx- 
imated as a constant density. Further work has confirmed 
the basic picture (Afshordi & ToUey 2008; McDonald 2008; 
Taruya et al. 2008). 

Slosar et al. (2008) have applied eq. (2) to constrain 
the value of /nl using a compilation of large-scale struc- 
ture data. They find -29 < /nl < +69 (at 95% confi- 
dence level). These limits are competitive with those from 
WMAP5, -9 < /nl < +111 (Komatsu et al. 2008) and 
-4 < /nl < 80 (Smith et al. 2009), demonstrating the 
promise of the method. Future all sky surveys could achieve 
constraints of the order of /nl ~ 5—10 (Dalai et al. 2008; Af- 
shordi & ToUey 2008; McDonald 2008; Carbone et al. 2008), 
assuming one knows how to extract maximum information 
from the data (see, e.g., Slosar 2008). In fact, with sufficient 
high density of tracers it should be possible to circumvent 
the sampling variance (which is a serious issue since the 
non-Gaussian effect is strongest on the largest scales) and 
alleviate degeneracies with other cosmological parameters, 
thereby allowing for a potentially huge gains (Seljak 2008). 

Still, in order to fully exploit the potential of forth- 



coming large-scale surveys, the method needs to be tested 
with large numerical simulations. Thus far, eq. (2) has been 
validated only using the halo-matter cross-power spectrum 
(Dalai et al. 2008) and only on very large scales, so its accu- 
racy remains uncertain. It is important to measure the effect 
in the auto-correlation of dark matter haloes, since the lat- 
ter gives the strongest constraint on /nl (Slosar et al. 2008). 
It is also important to extend the analysis to smaller scales, 
where the peak-background split breaks down, as well as to 
loss biased haloes. The purpose of this paper is to address 
these issues in more detail. We begin with a brief descrip- 
tion of the N-body simulations against which we calibrate 
the theory (§2). Next, wc discuss the mass function and bias 
of the corresponding halo catalogues and demonstrate the 
importance of including a scale-independent bias correction 
in the comparison with the simulations (§3). The main body 
of the paper is §4, where we study in detail the impact of lo- 
cal non-Gaussianity on the halo-matter and halo-halo power 
spectrum. We conclude with a discussion of the results in §5. 



2 THE N-BODY SIMULATIONS 

Investigating the scale-dependence of the halo bias requires 
simulations large enough so that many long wavelength 
modes are sampled. At the same time, the simulations 
should resolve dark matter haloes hosting luminous red 
galaxies (LRGs) or quasars (QSOs), so that one can con- 
struct halo samples whose statistical properties mimic as 
closely as possible those of the real data. 

In this work, we use a series of large N-body simula- 
tions of the ACDM cosmology seeded with Gaussian and 
non-Gaussian initial conditions. The non-Gaussianity is of 
the "local" form, $ = -|- /nl(</'^ - (0^)), where $(a:) is 
the Bardeen potential. It is important to note that this lo- 
cal transformation is performed before irmltiplication by the 
matter transfer function. T{k) is computed with CMBFAST 
(Seljak & Zaldarriaga 1996) for the WMAP5 best-fitting pa- 
rameters (Komatsu et al. 2008) : h = 0.7, f2m = 0.279, 
Qh = 0.0462, Us = 0.96 and a normalisation of the curva- 
ture perturbations = 2.21 x 10"® (at k = 0.02Mpc~^) 
which gives as « 0.81. Five sets of three 1024^ simulations, 
each of which has /nl = 0, ±100, were run with the N-body 
code GADGET2 (Springel 2005). We used the same Gaus- 
sian random seed field (j) in each set of runs so as to minimise 
the sampling variance. We also explored lower values of /nl 
and ran 2 realisations for each of the non-Gaussian models 
characterized by /nl = ±30 and ±10. In all cases the box 
size is 1600 /i~^Mpc with a force resolution of 0.04 times the 
mean interparticle distance. The particle mass of these sim- 
ulations thus is 3.0 X 10^^ Mq/H, enough to resolve haloes 
down to 10" Mq/H. 

Haloes were identified using the MPI parallelised ver- 
sion of the AHF halo finder which is based on the spherical 
overdensity (SO) finder developed by Gill, Knebe & Gibson 
(2004). AHF estimates the local density around each halo 
centre using a top-hat aperture. The virial mass M is de- 
fined by the radius at which the inner overdensity exceeds 
Avir(z) times the background density p{z). Note that Avir(«) 
is an increasing function of redshift (Avir « 340 at « = 0). 
We discard poorly resolved haloes and only study those con- 
taining at least 34 particles to reduce the error in the mass 
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Figure 1. Top panel : multiplicity function f{u, 0) for the Gaus- 
sian simulations. Different symbols refer to different redshifts as 
indicated. Results are shown relative to the Sheth-Tormen fitting 
formula to emphasise deviation from the latter. Bottom panel : 
ratio between the non-Gaussian and the fiducial Gaussian mass 
functions. The dotted and dotted-dashed curves are the theoret- 
ical prediction at 2 = and 2, which is based on an Edgeworth 
expansion of the dark matter probability distribution function 
(see text). In both panels, error bars denote Poisson errors. For 
illustration, M = 10^^ Mq/H corresponds to i' = 3.2, 5.2, 7.7 at 
redshift z = 0, 1 and 2, respectively. Similarly, M = 10^'* Mq/H 
and 10^3 Mq/H correspond to u = 1.9, 3, 4.5 and 1.2, 1.9, 2.9 
respectively. 



estimate (Warren et al. 2006). This implies a lower mass 
limit M = 10^'' M0//1 which is about the typical mass of 
QSO-hosting haloes at 1 < 2 < 2 (e.g. Porciani & Norberg 
2006) and is a few times smaller than the mass of haloes 
harbouring LRGs in SDSS (Mandelbaum etal. 2006). 



3 HALO MASS FUNCTION AND BIAS 
3.1 Multiplicity function 

Analytic arguments based on the Press-Schechter theory 
(Press-Schechter 1974; Bond et al. 1991; Sheth & Tormen 
1999) predict that the halo mass function n(M, z) is en- 
tirely specified by the distribution yf{v) of first-crossings, 
or multiplicity function 



M 



2 n{M, 2) dlnM 
p d In 



(3) 



The peak height v{M,z) = Sc{z)/cr{M), where 5c{z) « 
1.68D(0)/D(2) is the critical linear overdensity for collapse 
(assumed spherical throughout this paper), is the typical 
amplitude of fiuctuations that produce haloes of mass M 
by redshift 2. A characteristic mass for clustering, M*(2), 
can then be defined through v{M,z) — 1. For the present 
cosmology, M*(0) « 3.5 x 10^^ Mq/Zi. 

The top panel of Fig. 1 shows the multiplicity function 
of the SO haloes extracted from the Gaussian simulations 



at redshift 2 = 0.3, 0.5, 1 and 2. The numerical data are 
plotted with respect to the Sheth-Tormen function (Sheth- 
Tormen 1999) to emphasise the large deviation from the lat- 
ter. This departure is, however, not really surprising since 
the Sheth-Tormen formula is a fit to the mass function of 
Friends-of-Friends haloes (extracted from the GIF simula- 
tions, see Kauffmann et al. 1999). We have not attempted 
to fit the multiplicity function of our SO haloes given the 
limited volume and dynamic range of our simulations. In- 
stead, we have found more useful to assess whether the im- 
pact of local non-Gaussianity on the halo mass function is 
consistent with theoretical expectations. 

To test this we have plotted the ratio /(;/, fNh)/ fiy, 0) 
in the bottom panel of Fig. 1 for the simulations with 
/nl = ±100. The presence of primordial non-Gaussianity 
enhances or suppresses the high peak tail of the multiplic- 
ity function depending on the sign of /nl- As recognised 
in previous papers (e.g. Matarrese et al. 2000; Sefusatti et 
al. 2007), despite the lack of a reliable Gaussian mass func- 
tion, deviations from Gaussianity can be modelled analyt- 
ically using the Press-Schechter formalism. Here we follow 
the simple extension introduced by LoVerde et al. (2008; 
see also Chiu et al. 1998) and replace the Gaussian proba- 
bility distribution function (PDF) of the density field by the 
generic Edgeworth expansion (e.g., Scherrer & Bertschinger 
1991; Juskiewicz et al. 1995). Neglecting cumulants other 
than the skewness S'i{M) — {S'm) / (S^i)"^ and truncating the 
series expansion at Sa, the non- Gaussian correction factor 
reads (LoVerde et al. 2008) 

/('^../nl) 



/(^,0) 



1 + 



1.3 X 1 d(ag3) / 1\ 



(4) 



after integration over regions above the critical density 
for collapse. Note that we have omitted the explicit red- 
shift dependence. Strictly speaking however, the ratio 
/('^i /nl)//(j^, 0) depends distinctly upon the variables M 
(or v) and 2 due to the presence of (tS3,{M). Our notation 
is motivated by the fact that the measured non-Gaussian 
correction, as plotted in the bottom panel of Fig. 1, appears 
to depend mostly on the peak height. 

Equation (4) requires knowledge of the skewness Ss,{M) 
of the smoothed density field 5m, which we compute analyt- 
ically using the relation (see Appendix §A) 



a*S-i{M) = 



(27r2 



dfci kia{M,kx)P4,{k-i) 



(5) 







X / dk2kia{M,k2)P4,{k2) 
'0 

+1 

X / dfia(M,k) 1 + 2 



P^{k2) 



where k'^ = kf + k2 + 2/ifcifc2, Pcf,{k) is the power spectrum 
of linear curvature perturbations in the matter-dominated 
era. 



q(M, k) = 



30, 



772 



D{z)k^T{k)W{M,k) 



(6) 



and W{M,k) is a (spherically symmetric) window function 
of characteristic mass scale M . Over the mass range probed 
by our simulations, lO" <^M ^5 x 10^*^ M0//1, aSsiM) is 
a monotonic decreasing function of Al that varies in the 
narrow range 3 — 3.3 x 10~^/nl for the top-hat filter 
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Figure 2. Halo bias as a function of wavenumber. Results are 
shown at redshift z = 0.3, 0.5, 1, 1.4 and 2 (from bottom to top) 
for haloes with mass above 2 X lO'^^ Mq/Zi. Filled and empty 
symbols represent the bias estimators by^^^ = -y/ Phh/ Pmm and 
f'mh = ^mh/^mm, respectively. The bins are equally spaced in 
logarithmic space with a bin width A log A: = 0.1. Measurements 
of fehh have been slightly shifted horizontally for clarity. The hor- 
izontal lines indicate our estimate of the linear bias b{M) (see 
text). Notice that Phh(^) is corrected for shot-noise. 



assumed here. Furthermore, the aSz term dominates the 
total contribution to the non-Gaussian correction when the 
peak height is v ^2. 

The resulting non-Gaussian correction is plotted in the 
bottom panel of Fig. 1 for two different redshifts, z = Q (dot- 
ted) and 2 (dotted-dashed). The truncated expansion eq. (4) 
agrees reasonably well with the numerical data, suggesting 
thereby that cumulants higher than S3 may not be impor- 
tant in the range of mass and redshift considered here. Note 
also that for positive /nl the mass function is enhanced 
more at the high mass end and that this is similar to an 
increase in the amplitude of fluctuations ag. Hence, /nl is 
somewhat degenerate with erg since, in both cases, the ef- 
fect increases with mass (compare with Fig. 3 of Mandel- 
baum & Seljak 2007 for instance). However, &t u = 3.2 (i.e. 
M = 10^-^Mq//i at z = 0) the increase in mass function 
for /nl = 100 is 15%, which corresponds to less than 0.01 
change in erg. Therefore, given the current uncertainties in 
the cluster abundance (which translate into 0.03 error on 
ag, Vikhlinin et al. 2008), the prospects of using halo mass 
function to place competitive limits on /nl with the current 
data are small. 



3.2 Linear bias 

Having checked that the level of non-Gaussianity in the mass 
function is consistent with simple theoretical expectations, 
we now turn to the clustering of dark matter haloes. 

We interpolate the dark matter particles and halo cen- 
tres onto a regular cubical mesh. The resulting dark matter 
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Figure 3. Top panel : Non-Gaussian correction /3m(fc,/NL) = 
APmm(fc, /NL)/^'mm(fc, 0) to the matter power spectrum that 
originates from primordial non-Gaussianity of the local type. Re- 
sults are shown at redshift z = and 2 for /nl = ±100. The 
dashed curves indicate the prediction from a leading-order per- 
turbative expansion. Bottom panel : Non-Gaussian bias correction 
for the haloes of mass M > 2 X 10^^ Mq/H extracted from the 
snapshot at 2; = 0.5 (filled symbols). The solid curve represents 
our theoretical model eq. (9). The dashed, dotted and dashed- 
dotted curves show the three separate contributions that arise at 
first order in /nl- Our theoretical scaling agrees very well with 
the data for k ^0.05 /iMpc~l. 



and halo fluctuation flelds, (5m (fc) and Sh{k), are then Fourier 
transformed to yield the matter-matter, halo-matter and 
halo-halo power spectra Pmm(fc), -Pmh(fc) and Phh(fc), respec- 
tively. Notice that the power spectra are computed on a 512'' 
grid to reduce the computational expenses. Still, the Nyquist 
wavenumber is sufficiently large, ~ 1 /iMpc~^, to allow for 
an accurate measurement of the power in wavemodes of am- 
plitude k <0.1 /iMpc~^. As we will see shortly, the impact 
of local non-Gaussianity is negligible at fc = 0.1 /iMpc~^, 
but increases rapidly with decreasing wavenumber. 

On linear scales, the halo bias b{k) — 5h(fc)/5m(fc) ap- 
proaches a constant, albeit mass- dependent value b{M). The 
linear halo bias b{M) needs to be measured accurately as it 
controls the strength of the scale-dependent bias correction 
induced by local non-Gaussianity. To proceed, we may con- 
sider the following estimates of fe(fc) for a given halo sample. 



&hh(fc) 



Phh(fc) 

(k)' 



bmh(fc) 



Pmh(fc) 

(k) 



(7) 



In the following, we will always correct the halo power spec- 
trum for shot-noise, which we assume to be l/fih if dark mat- 
ter haloes are a Poisson sampling of some continuous field. 
While this discreteness correction is negligible for Pmm(fc) 
and Pmh(fc) due to the large number of dark matter parti- 
cles, it can be quite significant for Phh(fc). 

In Fig. 2, the result of measuring 6mh(fc) and 6hh(fc) in 
the Gaussian simulations is shown at various redshifts for the 
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haloes of mass M > 2 x 10^"^ Mq/Zi. Error bars indicate the 
scatter among the various realisations. Except for the most 
biased sample, &mh(fc) and bhh{k) are nearly constant and 
agree well with each other when the wavenumber varies in 
the "linear" range ~ 0.005 — 0.05 /iMpc~^. On these scales, 
the slight offset between bmh{k) and bhh{k) suggests that the 
shot-noise correction 1/nh might be too large for the low bias 
haloes and too small for the highest bias halo. It is worth 
pointing out that the hypothesis of shot noise being l/fih 
for the dark matter haloes remains unproven (McDonald 
2008) and it is an issue worth exploring further. Here we 
will be mostly looking at ratios of power spectra with and 
without non-Gaussianity, so this is less of an issue. Both bias 
quantities feature some scale-dependence on smaller scales, 
k Jj0.05 — 0.1 /iMpc^^. This is best seen in the most biased 
sample. We will use 6mh(fc) as a proxy for the linear halo bias 
since it is less sensitive to shot-noise. In Fig. 2, the horizontal 
lines indicate our fit to b{M) obtained from the measurement 
of &mh(fc) at wavenumber 0.005 < fc < 0.05 /iMpc"^. 

3.3 Non-gaussian bias shift 

As shown in Dalai et al. (2008), Matarrese & Verde (2008) 
and Slosar et al. (2008), local non-Gaussianity gives rise 
to the scale-dependent bias correction eq. (2). However, at 
the lowest order there are two additional, albeit relatively 
smaller, corrections which arise from the dependence of both 
the halo number density n{M, z) and the matter power spec- 
trum Pmm on /nl. As we will see shortly, the inclusion of 
these extra terms substantially improves the comparison be- 
tween the theory and the simulations. 

Firstly, assuming the peak-background split holds, the 
change in the mean number density of haloes induces a scale- 
independent shift which we denote by Abi(/NL). The exis- 
tence of such a term was noted in Slosar et al. (2008) and 
Afshordi & ToUey (2008). Using the non-Gaussian fractional 
correction eq. (4) (which is not universal), this contribution 
reads 



10 F 



(8) 



d m 



dVSa) 
din 



This approximation should work reasonably well for moder- 
ate values of the peak height, u <^4, for which the formula of 
Lo Verde et al. (2008) matches well our data (see Fig. 1). It 
is worth noticing that A6i(/nl) has a sign opposite to that 
of /nl (because the bias decreases when the mass function 
goes up). In practice, to estimate A6i(/nl) for a given halo 
sample, we evaluate a S3 and v at the scale corresponding to 
the average halo mass M of the sample. Furthermore, since 
we consider only first order corrections to the Gaussian bias, 
we set f{v,Q) = f{y,fNh) in the above expression so that 
Afei is truly first order in /nl. 

Secondly, primordial non-Gaussianity affects the mat- 
ter power spectrum as positive values of /nl tend to in- 
crease the small-scale power (Scoccimarro et al. 2004; Grossi 
et al. 2008; Taruya et al. 2008). For /nl ~ O(IO^), the 
magnitude of this correction is at a per cent level in the 
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Figure 4. Top panel : A comparison between the auto-power 
spectrum with and without the shot-noise correction. PhhC^, /nl) 
is measured at 2 = 1 for haloes of mass 

Prom top to bottom, the various symbols represent the simula- 
tion results with /nl = +100 (blue), (green) and -100 (red). 
The linear bias of this sample is b{M) ^ 2.5. Bottom panel : 
Phh{k, /nl)/ Phh{k, 0) as a function of wavenumber. The dashed 
curves denote the theoretical prediction (see text). In both pan- 
els, measurements without the shot-noise correction have been 
shifted horizontally for clarity. 



weakly nonlinear regime k <^0.1 /iMpc^^. In order to illus- 
trate this effect, the top panel of Fig. 3 displays the de- 
viation /3m(fc,/NL) = APnim(fc, /NL)/f'mm(fc, /nl = 0) that 

arises from the presence of primordial non-Gaussianity of the 
local type. The symbols show the result of measuring this 
ratio from the snapshots at redshift z = and 2, whereas the 
dashed curves show the prediction from one-loop perturba- 
tion theory (Taruya et al. 2008; see also Appendix §A). As 
we can see, leading order perturbation theory (FT) provides 
an excellent description of the effect over the wavenumbers 
of interest, k <^0.1 /iMpc~^. At z — 0, one-loop FT overes- 
timates the non-Gaussian correction by ~ 15 per cent for 
k = 0.1 h~^Mpc and it is possible the agreement could 
be improved further using renormalised perturbation the- 
ory (see, e.g., Crocce & Scoccimarro 2008). 

Summarizing, local non-Gaussianity adds a correction 
Afe(fc, /NL) to the bias b{k) of dark matter haloes that can 
be written as 



A6(fc,/N 



Afe«(fc, /nl) + A6i(/nl) + b{M)l3m{k, /nl) (9) 



at first order in /nl- The bottom panel of Fig. 3 illustrates 
the relative contribution of these terms for haloes of mass 
M > 2 X 10" Mq/H identified at redshift z = 0.5. The 
solid curve shows the total non-Gaussian bias A6(fe,/NL). 
Considering only the scale-dependent shift Afe„ leads to 
an apparent suppression of the effect in simulations rela- 
tive to the theory. Including the scale-independent correc- 
tion A6i considerably improves the agreement at wavenum- 
bers k ;^0.05 /iMpc~^. Finally, adding the scale-dependent 
term b{M)f3m further adjusts the match at small scale 
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k >0.05 /iMpc ^ by making the non-Gaussian bias shift less 
negative. 



4 RESULTS 

In order to quantify the effect of non-Gaussianity on the 
halo bias, we will consider the ratios ^ 



(10) 



Pmh(fc,/NL) ^ Ab(fc, /nl) 

Pmh(A;,0) 6(M) 

J^h(fc,/NL) ^ / A&(fc,/NL] 

Phh(fc,0) ^ 6(M) 

Moreover, we shall also quantify the departure from the 
theory as a function of wavemode amplitude with the ra- 
tio Ab^/A6*. Here, A6* is the non-Gaussian bias correction 
measured from the simulation wherecis A6* is the theoretical 
scaling eq. (9). 

Before proceeding we look at the effect of the shot- 
noise correction on the measurement of the non-Gaussian 
bias A6(fc,/NL). In Fig. 4, the averaged halo power spec- 
trum and the ratio Phh(fe, 0)/Phh(A:, /nl) are shown before 
and after applying the discreteness correction. Error bars 
represent the scatter among the realisations. The bias and 
the number density of the halo sample considered here is 
b{M) « 2.5 and fih « 10"" ft^Mpc"^, respectively. 

As we can see, the shot-noise can have a non-negligible 
effect on the largest scales, specially for the haloes ex- 
tracted from the simulations with /nl = —100 for which 
the large scale power crosses zero on very large scales. For 
this particular sample, the shot-noise correction enhances 
the measurement of A6(A;, /nl) by 10-15 per cent at scales 
k ^0.03 feMpc"^, regardless of the sign of /nl- While the 
exact amount of correction depends upon the bias and the 
number density of the halo sample under consideration, it 
is clear that any attempt to measure A6(fc, /nl) at the few 
per cent level must include the discreteness correction. 

4.1 Non-Gaussian bicis from the halo- halo and 
halo-matter power spectra 

We have measured power spectra for a range of halo masses 
and redshifts, covering the relevant range of statistical prop- 
erties corresponding to the available data sets of galaxies or 
quasar populations with different luminosities and bias. The 
results are summarised in Figures 5 and 6, where the aver- 
aged A6/6 and (1 -I- Ab/b)'^ — 1 arc plotted as a function of 
wavenumber. The deviation from the theoretical prediction, 
Afe^/Afe*, is also shown at the bottom of each panel. The 
shaded region indicates a deviation less than 20 per cent. 
To reduce the impact of sampling variance, we first compute 
the ratios Pmh{k, /nl) / Pinh{k,Q) and Phhik, /nl) / Phh{k,Q) 
for each realisation, and then average over the realisations 
(see, e.g.. Smith et al. 2007). We note that reversing the se- 
quence of operations, i.e. taking the ratio of averaged power 
spectra, gives very similar average values. Error bars denote 



^ Strictly speaking, Phh{k, fNh)/Phh{k,0) is equal to [1 -|- (A^ -l- 
Ai)/6]2 + _ 1, which differs from eq. (10) by /3m -|- 0{l3^). In 
what follows however, we will use (1 -|- A6/6)^ — 1 for notational 
convenience. 



the scatter around the mean and, therefore, may underesti- 
mate the true errors since they are computed from a small 
number of realisations. 

As we can see, the theoretical prediction provides a very 
good description of the simulations at small wavenumber 
k <,0.0S /iMpc"\ but the ratio Ab"/Ab* differs significantly 
from unity at larger wavenumbers. The exact amount of de- 
viation depends weakly on the sign of /nl. For moderately 
biased haloes with 2 < b{M) < 3 the theory approaches 
the numerical results already on scale k ;^0.05 /iMpc~^. For 
the highly biased samples 6 > 3, the theory overpredicts 
the effect seen in simulations on all scales, but somewhat 
more on smaller scales, although in the high bias limit the 
numerical data is noisier due to the very low number den- 
sity of haloes. It is worth noticing that, at the largest scales 
k ^0.005 hMpc~^, the cross-power spectrum Pmh(A;,/NL = 
—100) goes negative while Phh{k, /nl = —100) remains posi- 
tive and even increases, in good agreement with the analytic 
prediction. 

We suspect these deviations at high wavenumber are 
mostly due to the breakdown of the peak-background split 
approximation which was used in the derivation of the scale- 
dependent bias term Ab^ in Slosar et al. (2008). For this 
approximation to be valid one assumes that the long wave- 
length modes act as a homogeneous change of the back- 
ground, from which the effect of the non-Gaussianity is com- 
puted by comparing it to the local rescaling of the fluctua- 
tion amplitude. Clearly this assumption breaks down once 
the wavelength of the mode becomes small. Uncertainties 
in the scale- independent correction also affect Ab" /Ab*. In 
this paper we use analytic predictions based on equation 
(8), but we could also treat the scale- independent bias as a 
free parameter that we fit to the data, as done in the ac- 
tual data analysis of Slosar et al. 2008. For example, a ~ 20 
per cent smaller (larger) A6i at 6(M) ^3 (6(M) ^3) would 
noticeably improve the convergence at large k. Finally, no- 
tice that the auto- and cross-power spectrum of haloes give 
comparable results at all but the (poorly sampled) largest 
scales, where sampling variance prevents us from making 
any conclusions. This confirms the validity of the analysis 
in Slosar et al. (2008), where this effect was applied to the 
auto-correlations of galaxies and quasars. 

To assess the extent to which the agreement between 
simulation and theory depends upon the halo mass and bias. 
Figures 6 and 7 further explore the effect in the low and high 
redshift outputs. In Fig. 6, the non-Gaussian bias is shown 
for haloes that correspond more closely to the quasars used 
by Slosar et al. (2008), which are at 2 = 1.8 and with b = 2.7. 
Our halo samples span a similar redshift range, 1.4 < z < 2. 
However, the mass cut 10^^ < M < 2 x 10^^ Mq/Ii gives 
larger values of the bias, 3 <^b{M) <^5, suggesting that the 
quasars are hosted by haloes (slightly) less massive than 
10^^ Mq//i (unresolved in our simulations). As can be seen, 
the correction factor Ab^/Ab* is similar to that of the sam- 
ples at high redshift z > 1 (cf. Fig. 5). 

For the redshift outputs z < 0.5, the relatively large 
number of dark matter haloes allows us to split the cat- 
alogues into several non-overlapping subsamples having a 
number density Uh — W~* h'^Mpc'^ . For these snapshots, 
we consider the mass bins 10^^ < M < 1.6 x 10^^ Mq/H, 
1.6 X 10" < M <3 X 10" Mo/Zi and M > 3 x lO" Mq/H. 
Results are shown in Fig. 6. An increase in the ratio Ab'^/Ab* 
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Figure 5. Non-Gaussian bias correction measured in the simulations at various redshifts for haloes of mass 

(colors as in Fig. 4). In each panel, the upper plot shows the ratio Phh{k, fNh)/ Phhi^^yO) — 1 (dotted curves, empty symbols) and 
^mhCi^i fNh)/ Pmhi^y 0) — 1 (solid curves, filled symbols). The error bars represent the scatter among 5 realisations. The respective output 
redshift and linear halo bias are also quoted. The bottom of each panel displays the departure from the theoretical prediction, Afe*/A6* 
(see text). The shaded area indicates the domain where the deviation is less than 20 percent. The theory agrees reasonably well with the 
measurements at wavenumber k ^0.03 /tMpc^^. 
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Figure 6. Same as Fig. 5 but for haloes extracted from tlie simulation outputs at z = 1.4, 1.7 and 2, with a mass in the range 
1 < M < 2 X 10^3. The halo sample at z = 1.4 is close to the QSO sample used by Slosar et al. (2008), for which z = 1.8 and b = 2.7. 



as a function of wavenumber followed by a change of sign 
can also be seen in these low biased samples in spite of the 
noisier data. Notice that the linear halo bias is in the range 
1 <p{M) <.2. In particular, the z = haloes with mass 
10^^ < M < 1.6 X lO" Mq//i constitute an almost un- 
biased sample of the density field, with b{M) ~ 1.10. At 
scales k <0.02 /iMpc~^ there is some evidence that the non- 
Gaussian bias correction measured in the low biased samples 
may be larger than the theoretical expectation. Still, the bias 
shift is quite small for b{M) = 1.10, in agreement with the 
theoretical prediction that the effect vanishes for b{M) = 1 
assuming the Eulerian bias prescription b{M) = 1-1- bi^{Ad) 
(where is the Lagrangian bias) used in eq. (9). Unfor- 

tunately, our simulations do not have sufficient mass resolu- 
tion to resolve anti-biased haloes with b{M) < 1, for which 
theoretical predictions based on the peak-background split 
suggest the sign of the scale-dependent contribution A6k is 
reversed. 

We have not examined the behaviour of Ab^/Ab' at 
k > 0.1 /iMpc~^ since the effect is already quite small there 



and nonlinear bias due to galaxy evolution effects domi- 
nates. Most of the information on the non-Gaussian bias 
comes from measurements at large scale k ;^0.03 /iMpc~^ 
(see Slosar et al. 2008), where the theoretical model and the 
numerical data agree reasonably well with each other over 
the relevant range 2 < 6 < 3. 

To reduce the scatter in the measurement of Ab{k, /nl), 
we can increase the bin width A log k so as to increase 
the number of independent modes. In Fig. 8, the ratio 
Ab^/Ab' is shown as a function of the linear halo bias for 
three equally spaced logarithmic interval spanning the range 
0.0045 < fc < 0.035 /iMpc"^ (e.g. Alogfe = 0.3). The data 
points are harvested from several outputs spanning the red- 
shift range < 2; < 2 (i.e., the snapshot redshifts are z = 0, 
0.3, 0.5, 1, 1.4, 1.7 and 2). The squares and triangles repre- 
sent the deviation from the theoretical prediction obtained 
by taking ratios of Pmh and Phh, respectively. Filled sym- 
bols show results for the non-Gaussian simulations with 
/nl ~ +100. The error bars indicate our jackknife error 
estimates on the average Ab^/Ab'. 
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Figure 7. Same as Fig. 5 but for haloes extracted from the simulation outputs at z = and 0.3, with a mass in the range 1 < M < 
1.6 X 10^'' (upper panels) and 1.6 < M < 3 X 10^'' Mq/Zi (lower panels). The sample with z = 0.28 and b = 1.65 roughly corresponds to 
LRG sample used by Slosar et al. (2008). 



The non-Gaussian bias shift of the low biased samples, 
& < 2, appears to deviate from the theory. Equation (9) may 
thus need correction when the linear bias gets lower than 
<1.5 — 2. The scale-dependent bias is quite large around 
b ~ 1.3, although for b ~ 1 the effect does appear to vanish 
on the largest scales as expected (see the upper left panel of 
Fig. 7). For the haloes with b(M) >2, the theory matches 
the non-Gaussian bias correction for k ^0.01 /iMpc~^. One 
would need even larger simulation boxes than used here to 
properly sample the largest scales. For k ^0.01 /iMpc~^, 
the effect is slightly suppressed compared to the theoretical 
prediction, but this may plausibly arise from uncertainties 
in the magnitude of the theoretical scale-independent shift 
A&i. Indeed, we have found that a 20 per cent increase in Abi 
considerably improves the agreement for k ;^0.05 /iMpc~^ 
and, at the same time, is still consistent with the measured 
fractional change in the multiplicity function (see Fig. 1). 
Finally, note that that haloes with similar bias also have 
a comparable scale-dependent bias due to non-gaussianity 



regardless of redshift. Hence, there is no need to introduce 
a second parameter such as redshift for the purpose of de- 
scribing these results. 

The only previous work with simulations along these 
lines is that of Dalai et al. (2008). These authors do not in- 
clude the scale-independent shift A6i nor the weaker correc- 
tion b{M)l3m induced by the matter power spectrum. Hence, 
Fig. 8 of their paper indeed shows l\b" / Ab^. This ratio ap- 
pears to increase with wavenumber (even though their data 
points do not extend beyond 0.03 ZiMpc"^), while the bot- 
tom panel of our Fig. 4 shows that Afe^/AbK is suppressed 
at high wavenumbers. Note, however, that their theoretical 
scale-dependent correction AbK(fe, /nl) does not include the 
matter transfer function. We found that, if the transfer func- 
tion were removed from eq. (2), Ab^/AbK would be enhanced 
rather than suppressed as one goes to higher wavenumber, 
in qualitative agreement with their findings. Since the non- 
Gaussianity is imprinted in the initial conditions prior to 
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Figure 8. Ratio of simulations to theoretical predictions, eq. (9), 
as a function of linear halo bias. The scale-dependent correction 
is calculated using all the wavemodes in the wavenumber ranges 
quoted on the figure. Squares and triangles indicate the value of 
the ratio Aft'/Ab* calculated from P^h a-nd Pmhi respectively. The 
measurements from the non-Gaussian simulations with /nl = 
-1-100 are marked as filled symbols. The error bars are computed 
from a jackknife estimate. 



Figure 9. Sensitivity of the non-Gaussian shift to the strength 
of the nonlinear parameter /nl. The ratio Ab{k, f^j^)/ Ab{k, /^l) 
is plotted as a function of wavenumber for various values of f^^^ 
and f^j^ spanning the range [-100,-1-100]. Symbols show results 
ai z = 0.3 (triangle) 0.5 (square), 1 (circle) and 2 (cross) for two 
different halo mass cuts : M > 2 X 10^'^ Mq/Zi (upper panels) 
and 1 < M < 2 X 10^^ Mq/H. The horizontal lines indicate the 
linear scaling /nl//nl- 



the evolution through matter and radiation domination, the 
transfer function must be included in the analysis. 

4.2 Scaling with /nl 

The quadratic term /nl0^ also induces second and higher 
order corrections to the effective bias shift A6(fc,/NL) 
which may become important at high wavenumber. To 
test for these high order terms, we explore in Fig. 9 the 
scaling of the non-Gaussian bias shift with the strength 
of the nonlinear parameter /nl. Symbols show the ratio 
A6(fc, /iJ/A^fc, /Sl) (which we abridge AfcC/iJ/AK/^J 
for shorthand convenience) as a function of wavenumber and 
redshift for several values of /^l and /^l spanning the range 
[-100, -f 100], as indicated in the figure. Note that the data 
points are obtained by averaging over two realisations only. 
The horizontal line indicates the value /nl//nl that should 
be reached if the non-Gaussian bias shift is linear in /nl. 

As we can see, there is less scatter in 
A6(±100)/A&(±30) than in A6(±30)/A6(±10) but, in 
both cases, the results are broadly consistent with the 
linear expectation /nl//nl- Furthermore, there is no 
significant dependence on the wavenumber, redshift or the 
halo mass cut. We conclude that the sensitivity of large 
scale structure bias should extend to smaller values of /nl 
as expected. 

4.3 Non-Gaussian bias in configuration space 

Thus far, we have investigated the impact of local non- 
Gaussianity on two-point statistics in Fourier space. It is 



also instructive to consider the two-point correlation ^(r) in 
configuration space, which is related to the power spectrum 
P{k) through 

1 

i{r) J dfcfcV(fc)io(fcr) , (11) 

where jo{x) is the zeroth spherical Bessel function. In prac- 
tice, since the simulation volume is a periodic cube, we com- 
pute the correlation from a discrete Fourier transform of the 
power spectrum. 

In Fig. 10, the result of measuring the auto and cross- 
correlation functions is shown at 0.3 < z < 1.5 for the 
mass cut M > 2 x 10^^ Mq/Zi. The width of the simulation 
box is large enough to sample wavemodes relevant to the 
baryon acoustic oscillations (BAG). The interesting feature 
of Fig. 10 is the correlation between the BAG and the broad- 
band power, which shows up differently in the correlation 
function than in the power spectrum. Local non-Gaussianity 
adds broadband power and, therefore, modulates the ampli- 
tude of the BAG and the position of zero-crossing. 



5 DISCUSSION AND CONCLUSIONS 

The scale dependence of clustering of biased tracers of the 
density field has emerged as a powerful method to constrain 
the amount of primordial non-Gaussianity of the local type. 
In this paper, we have measured the non-Gaussian bias cor- 
rection A&(fc, /nl) in the clustering of dark matter haloes 
extracted from a suite of large N-body simulations. In con- 
trast to previous work, we focus both on the halo-halo and 
halo-matter power spectrum. While we confirm the basic 
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Figure 10. Effect of local non-Gaussianity on the auto and cross- 
correlation function of haloes and dark matter, ^hh(^) ^^nd ?mh('')- 
Results are shown as a function of comoving separation r for the 
samples with M > 2 x 10^"^ Mq /h. The respective values of output 
redshift and linear halo bias arc also quoted. 



effect reported in Dalai et al. (2008), we emphasize the im- 
portance of including a scale-independent term A6i and, to 
a lesser extent, a contribution induced by the matter power 
spectrum b{M)l3-m, to the scale-dependent shift Ab^ when 
comparing the theoretical scaling to numerical simulations. 
The inclusion of these two first order corrections significantly 
improves the agreement at wavenumber k <P-1 /iMpc~^. 

The original analysis in Dalai et al. (2008) only used 
cross-power spectra from simulations, while the data analy- 
sis in Slosar et al. (2008) used mostly auto-power analysis. 
The two do not have to agree with each other if the haloes 
and dark matter do not trace each other on large scales, i.e. 
if there is stochasticity. While models with Gaussian initial 
conditions predict there is little stochasticity on large scales 
(Seljak & Warren 2004), this has not been shown explicitly 
for models with non-Gaussianity. Hence, one of the main mo- 
tivations for this work was to extract the non-Gaussianity 
effect from the auto-correlations. Measurements of the non- 
Gaussian bias correction obtained with the halo-halo or the 
halo-matter power spectrum are in a good agreement with 
each other, indicating that non-Gaussianity does not induce 
stochasticity and the predicted scaling applies equally well 
for the auto- and cross-power spectrum. The issue of stochas- 
ticity in non-gaussian models will be explored further in a 
future publication. 

For biased haloes (6 ^1.5), our results indicate that the 
simulated non- Gaussian bias converges towards the theoret- 
ical prediction for k ^0.03 /iMpc~^. At smaller scales, the 
effect depends on scale-independent bias. If it is ignored then 
the amplitude of the effect is suppressed relative to theory. If 
we include scale independent bias using analytic calculation 
this suppression is much smaller and in some cases goes in 
the opposite direction. Moreover, one could argue that scale- 
independent bias cannot be identified from the data alone, 
so one should fit for it and include it in the overall bias, as 
was done in Slosar et al. 2008. In this case the agreement 



between theory and simulations is improved further. Still, 
there is some evidence that for very biased haloes, & > 3, 
the effect is suppressed relative to theory even on very large 
scales. 

For the halo samples with b{M) ^1.5 there is some ev- 
idence that the actual bias exceeds the theory on all scales. 
Therefore, the proposed eq. (9) does not appear to be uni- 
versal, so care must be exercised when applied to the actual 
large scale structure data. It would be useful to verify eq. (9) 
on dark matter haloes which are anti-biased {h(M) < 1) rel- 
ative to the matter distribution, to see if the sign of the effect 
is reversed. One, however, needs a very large volume and a 
very high mass resolution, which prevents us from verifying 
the predictions in this regime with the current simulations. 

On the observational side, Slosar et al. (2008) have al- 
ready applied the method to a sample of highly biased LRGs 
and QSOs, with mean bias b{M) ~ 1.8 and 2.7, respectively. 
It is interesting to inspect how those constraints change in 
light of our analysis. Our results suggest that for these values 
of halo bias theory and simulations are largely in agreement 
on relevant scales: their constraints arise mostly from the 
measurement of the quasar power spectrum with b ~ 2.7 
at the largest angular scales, k <^0.005 /iMpc~^ and from 
LRGs with 6 ~ 1.8 at k ^0.01 /iMpc"\ As we see from Fig. 
8, theoretical predictions are in very good agreement with 
the simulations for these values of bias and scales. Hence, 
we thus expect their limits remain unchanged. 

Finally, we note that we have not considered other ef- 
fects that may also modify the predictions, such as redshift 
space distortions and merger bias. The latter can signifi- 
cantly weaken the predicted scale-dependent bias (Slosar et 
al. 2008). We plan to investigate these effects with simula- 
tions in the future. 
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APPENDIX A: PERTURBATION THEORY 
WITH LOCAL NON-GAUSSIANITY 

Al Skewness parameter 

In /nl non-Gaussiariity, the Fourier mode of the curvature 
perturbation (after matter-radiation equality) is given by 



$(fc) = <t>{k) + /n 



(f>{q)ct)(k- q) , 



(Al) 



where is the unperturbed Gaussian field with power spec- 
trum P4,(k) cc k^'~*. The primordial bispectrum of curva- 
ture perturbations is 



B<i>(fci, fca, fca) = 2/nl [P0(fci)P0(fc2) + 2 perms ; 



(A2) 



Hence, the three-point correlation of the Fourier modes of 
the smoothed matter density field, (5M(fc) = a(M, /c)$(fc), 
reads 



(5M(fa)^M(fe)^M(fc3)) 



(27r)'^ aiQ203 Bi,{ki, k2, ks) 
xSuiki + ki + ks) , (A3) 



where Oi — a{M,ki) for shorthand convenience. Here, So is 
the Dirac delta and the transfer function a{M,k) is given 
by eq. (6). Note that we have omitted the explicit red- 
shift dependence of 5m and a for brevity. The (connected) 
three-point function of 5m in configuration space is the 
Fourier transform of {5M{ki)5M{k2)dMik3)). In particular, 
the third-moment of the smoothed density field is 



{5li) 



d-^ki 



(27r)3 J (27r)3 J (27r) 
2/nl 



■{5M{ki)5M(k2)5M{k3)} 



f d'ki f 

'J i^^rJ' 



d^k2 



1 + 



P^ka 



(27t)3 

P^k-i 



aia2a3 P4,{ki)P^{k2) 



P4,(ki) P4>{k2) 



(A4) 



We have used the momentum conservation implied by the 
Dirac delta, i.e. ks = —ki — k2, to obtain the second line. 
Eq. (5) follows after taking advantage of the invaxiance un- 
der the exchange of fei with k2 and integrating out some of 
the angular variables. 



A2 Matter power spectrum 

Following Taruya et al. (2008), we estimate the non- 
Gaussian correction to the matter power spectrum in the 
weakly nonlinear range, k ^0.1 /iMpc~^, using perturba- 
tion theory. At the first order, the matter power spectrum 
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can be expressed as 

Pmm(fc,/NL) = D' (z) (k) + [P^''^\k, z) + P^'''\k , z)] 

+P(i''(fe,a;/NL) . (A5) 

Here, D{z) is the growth factor, Ph{k) is the linear power 
spectrum of the density field, 

u3 



P^'^'>{k,z) = 



D\z) 



98(27r)2 



da; Pl i^^) 



(A6) 



^ 

/ AjJ, P^{k^/l+3p^^^2|lx) 



P^''\k,z) = 



D\z) 



3a; + 7/i — lO/i x 

1 + a;2 - 2^.x 

k'P^k) 



252(27r)' 



Jo 



Ax Pl (fcx) 



(A7) 



12 

a;2 

3 / 2 



158 + lOOx - 42a; 



+ -L(x^-l)^(7a;^ + 2)ln 



1 + a; 



1 - X 



are the standard one-loop contributions in the case of Gaus- 
sian initial conditions (e.g. Goroff et al. 1986; Makino et al. 
1992; Jain & Bertschinger 1994) and 



P^'^\k,z;hi.) = /n 



2k^ 



'7(27r)2 



a(0, fc) / da; a; a(0, fcx) 





7/i — 10/t a; 



a(0,g) 
(A8) 



+ a;^ — 2fj,x 

X [P^(fe)P,^(A;a;) 2 perms ] 

where = A;^(l + a;^ — 2/ia;), is the leading-order cor- 
rection due to local non-Gaussianity which arises from 
the non-zero primordial bispectrum of curvature pertur- 
bations. a{0,k) is the function eq.(6) with filtering kernel 
W{0,k) = 1. The particular redshift dependence of these 
power spectra follows from the assumption of growing-mode 
initial conditions. The relative contribution /3m(fc,/NL) = 
APmm(fe, /NL)/Pmm(A;, 0) of local non-Gaussiauity thus is 

P(i2)(fc,z;/NL) 



/3m(fc, /n 



D2(z)PL(fc) +P(22)(fc,z) +P(13)(fc,2) 



(A9) 



at leading order. Notice that this ratio scales as oc D{z), 
so the effect of local non-Gaussianity on the matter power 
spectrum is largest at low redshift. 
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